Universal supercritical current distributions in 3D random resistor networks 
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Inspired by percolation-induced properties of nano-particle composites [Nan et al, Annu. Rev. 
Mater. Res. 40, 131 (2010)], and in particular the robust electrical property gains at and above 
percolation threshold, we revisit the classical setting of percolation in random resistor networks. 
We simulate and analyze scaling behaviors of the percolation-induced distribution of currents in 
the network, first confirming the celebrated power-law distribution of small currents. We further 
identify a robust exponential distribution in the large current tail that persists above but close to 
threshold, whereas the power-law scaling in the small currents does not. 



Motivated by questions about macroscopic electronic 
properties in material science, we consider the classic 
problem of bond percolation [l|, 0] in a three-dimensional 
L x L x L cubic lattice and the state of currents con- 
ducting through the percolatingbonds. In the seminal 
work of Dc Arcangelis et al. [1,13 j a hierarchical lattice 
model is given for the percolating backbone of the net- 
work, which yields a log-normal current distribution in 
the network; this model was later generalized by Lin et al. 

These models capture many features of the network 
such as the multifractal behavior of the current distri- 
bution (i.e., an infinite hierarchy of exponents is needed 
to describe the moments of the distribution). However, 
they failed to predict the power law distribution of small 
currents (Straley @ and Duering et al. @, 0)- Even 
more importantly, the previous models do not appear to 
adequately capture the distribution of bond currents that 
dominate measurable behavior. 

Current distributions have fundamental importance in 
material science. Low moments of the current distribu- 
tion dictate physically measurable properties, ej* the 
second moment describes the bulk conductance [3|, [9| , so 
that measurable scaling behavior at or above percolation 
threshold (cf. 0) is inherited from the current distribu- 
tion. Another illustration is in the study of breakdown 
of random media This critical network property 

motivated studies on the size and location of the largest 
current in the network 15Hl7|. We note the maximum 



current does not inform properties of the large current 
distribution. Here we analyze the entire current distri- 
bution, reproducing previous results about small currents 
while revealing scaling behaviors of the large current tail 
of the distribution that controls bulk response above per- 
colation threshold. 

Each edge of the lattice takes conductance 1 with prob- 
ability p, (we call such a conducting edge a bond and p 
is the bond density) , and conductance with probability 
1 —p. We seek to understand the relationship between the 
electrical and topological properties of the resulting bond 
percolation network. In order to model an externally- 
driven bulk electrical response, we consider two perfectly 
conducting L x L plates to be present at opposite ends 



of the cube, representing the sink and source of current 
(in response to either an external voltage drop or current 
source). 

For each realization of the random resistor network 
model, a breadth-first search algorithm is applied to 
find the union of percolating clusters that connect the two 
plates. The 2-core of this union is identified, with forced 
inclusion of the two plates, because only these bonds 
can carry nonzero currents under the defined electrical 
problem. Near the percolation threshold, the number of 
relevant bonds in the system is reduced by almost 90% 
by this preprocessing. More importantly, we observe im- 
proved numerical precision by this a priori elimination of 
many of the zero-current bonds. The Kirchhoff law is 
then solved on these percolating 2-cores with a standard 
linear solver, giving the current on each bond. A statis- 
tical description of the nework properties is obtained by 
averaging over many realizations (typically 1000 realiza- 
tions for each bond density and system size). 

The bond percolation threshold p c for an infinite 3D 
cubic lattice is p c = 0.2488 [2oj | . For p < p c , no per- 
colating cluster forms (with measurable probability) and 
thus current does not pass through the network. For 
convenience, let f(i) be the probability density function 
(PDF) of the current magnitudes across the population 
of current-carrying bonds (that is, ignoring zero-current 
bonds where present). Let h{x) be the corresponding 
PDF (again, restricted to non-zero currents) of the log- 
arithmic current X = ln(7). The two distributions are 
related by h(x) = f(e x )e x . 

The logarithmic current distribution h(x) and current 
distribution /(i) near (p = 0.25) and above (p = 0.29) 
the threshold are shown in Figure Q] for a unit voltage 
source. First, the logarithmic transformation of current 
exposes the small current region; the left panel in Figure 
[1] shows that the small currents obey a power- law distri- 
bution up to apeak in the current distribution, agreeing 
with Straley @ and Duering et al. @,H1- Second, for 
relatively large currents (i.e., to the right of the peak of 
the h(x) distribution), FigurcQ] (right) suggests exponen- 
tial f(i) distributions. This general shape of the current 
distribution persists as bond density p increases near and 
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above threshold. The implications are two-fold: there is 
an exponential large current tail at and above percola- 
tion, whereas the small current tail disintegrates above 
the critical bond density. We are cautious to rule out fi- 
nite size effects since the p = 0.25 system size is compara- 
ble to the correlation length [l[ . Specifically, near thresh- 
old, the correlation length £ scales as £ oc (p — p c )~ 09 
which we will include in our study below. 

While the small-current behavior is well understood 
and the maximum current has been studied extensively, 
little is known about the large current tail of the dis- 
tribution even though the large currents dominate bulk 
properties. E.g., the second moment of the current dis- 
tribution determines the bulk conductance. We focus 
here on the large current tail, revealing a robust exponen- 
tial distribution above but close to percolation threshold. 
The empirical densities of the current for a unit voltage 
source near percolation threshold {p = 0.25) and above 
threshold (p = 0.29) for different system sizes are shown 
in Figure [2] and Figure [3] respectively. Despite slightly 
larger noise at bond density p = 0.25, the straight lines 
in both figures point to exponential tails of the current 
distributions, and the rate of the exponential decay in- 
creases with the system size L. In order to meaningfully 
capture an externally-imposed voltage drop in the ther- 
modynamic limit (L — > oo), and to better understand the 
effect of system size on the current distribution, we carry 
out a finite-size scaling analysis on the distributions. Let 
/i(i) be the probability density function (PDF) of the 
current at system size L for a unit voltage source; then 
by properly rescaling /l(i) with L, we aim to eliminate 
the effect of the system size: 

L- u f L {L- v i) = r(i), (i) 

where f°°(i) is a function independent of L, and u and v 
arc critical exponents to be determined. By tuning u and 
V, we confirm that the densities for different system sizes 
collapse onto a single curve with u = 1 and v = 1 (see 
the insets of Figure [5] and Figure [3]). This is the expected 
result for a material with bulk conductance: the total 
resistance per unit cross-sectional area increases ~ L. In 
other words, f°°(i) is the limiting current distribution 
for system size L and external voltage source V* = L 
(made dimcnsionless by the definition of bonds as unit 
conductances). Therefore in a finite system the current 
density for a unit voltage source scales as: 

f L (i)=Lf°°(Li). (2) 

Note that equation ([2} requires that any multifractal 
property of the overall current distribution comes from 
the small current tail, since the k th moment of the 
large currents described by this finite-size scaling is a 
simple scaling function of L. Specifically, 

/•OO />oo 

M k = i k f L (i)di= i k Lf°°{Li)di^L~ k . (3) 
Jo Jo 



To confirm this simple scaling form of the moments, we 
compute the first several sample moments at bond den- 
sity 0.25 for varying system sizes and plot them in Figure 
IH The sample moments are calculated as Mk = j? J2b *6> 
where N is the number of bonds with nonzero current 
and the sum of bond currents i& is taken over all current- 
carrying bonds, b. The scaling forms of the moments are 
not exactly the same as Equation J3]) due to the mul- 
tifractal property of small currents; however for large 
moments the exponential tail of the current distribution 
becomes dominant and thus the scaling relationship ap- 
proaches Equation (J3j> - 

To confirm the universal nature of the exponential tail 
in the large values of the current distribution, we now 
focus on a unit current source flowing between the two 
plates. Current distributions for a unit voltage source 
differ from those for a unit current source by a factor of 
the bulk conductance C, because of the linearity of the 
system. Formally, denoting the density of the current dis- 
tribution for the unit current source by <?l(«), the density 
fh{i) of the current distribution for a unit voltage source 
can be written as: = gL,(i/C)/C. Then assuming 

<?l(*) ~ e _Al for large currents yields 

f L {i) ~exp-^, (4) 

for large currents. Therefore, the rate of the exponen- 
tial tail of being independent of the bond density 
(over an observed range) indicates a linear relationship 
between the rate of decay of the exponential tail of gt{i) 
and the bulk conductance, and vice versa. 

It would seem surprising for these two properties to 
scale linearly with one another. To quantify this relation- 
ship, exponential distributions are fitted to the current 
distributions for a unit current source at various bond 
densities and their rates are plotted against the bond 
densities in Figure [5j along with the scaling behavior of 
the bulk conductance. These results indicate that the 
rate of the exponential decay A(p) and the bulk conduc- 
tance C (p) have similar scaling behaviors with respect to 
bond density near and above the percolation threshold, 
which can be deduced from the figure: 

A(p)oc(p- Pc ) 1 - 969±0 - 077 , (5) 
CWccb-ft) 1 - 91110 - 012 (6) 

Both the underlying exponential distribution and the 
power-law scalings become invalid further from the 
threshold as the system saturates. We notice their clear 
breakdown at p = 0.35 and beyond (as indicated in the 
supplemental materials plj]). 

Finally we superimpose all the limiting current distri- 
butions f°°(i) = L^ 1 /i(L _1 i) for a wide range of system 
sizes L and bond densities p in Figure El with colors rep- 
resenting the ratio L/{p — p c )~ 9 as an indicator of the 
extent to which our result is affected by the finite size ef- 
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X=ln(l) I 

FIG. 1. Distribution h(x) of the logarithmic current (left) and distribution f(i) of the current (right) near threshold (p — 0.25) 
and above threshold (p = 0.29) in 100 x 100 x 100 random resistor networks. 
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FIG. 2. (Color online) Empirical probability density function 
of the current in an L x L x L cubic lattice at bond density 
0.25 (which we note is quite close to the percolation threshold 
0.2488) for various system sizes L. A constant unit voltage 
is imposed across the system. The plot is an average of 1000 
realizations. The inset shows the same distributions rescaled 
by the system L with critical exponents u = 1 and v = 1. 
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FIG. 3. (Color online) Empirical probability density function 
of the current in an L x L x L cubic lattice at bond density 0.29 
(which is above the percolation threshold 0.2488) for various 
system sizes L. A constant unit voltage is imposed across the 
system. The plot is an average of 1000 realizations. The inset 
shows the same distributions rescaled by the system L with 
critical exponents u — 1 and v = 1. 



fcct. Despite the noise at low system-size-to-correlation- 
lcngth ratios, Figure H] demonstrates apparent conver- 
gence to a universal class of exponential distributions for 
the large currents both near and above threshold. The 
rate of the exponential tail does not depend on the bond 
density and the simple scaling form of the tail of the cur- 
rent distribution with the system size remains the same. 
Again this universality of current distributions only holds 
close to threshold. Far above threshold the current distri- 
bution approaches a delta function which is the current 



distribution at p = 1 where every bond on the paths per- 
pendicular to the two plates carries the same amount of 
current while all other bonds carry zero current. Again, 
we notice significant difference in the current distribu- 
tions starting at approximately p = 0.35 (as indicated in 
the supplemental materials |21[). 

Recall that the second moment of the current distribu- 
tion is related to the conductance of the network. Specifi- 
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FIG. 4. Scaling of the first 6 sample moments with respect to 
the system size L. The bond density is fixed at p = 0.25 and 
the system size L varies from 50 to 200. The fitted equations 
of the moments Mk are shown in the figure (cf. Equation ©). 
Curves are normalized by their values at L — 200. 
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FIG. 5. Scaling behaviors against bond density p — p c for, 
the bulk conductance (C(p)), the rate of the exponential tail 
of the current distribution (A(p)), and the number of bonds 
with current larger than 10 -3 (N(p)). The system size is fixed 
at L — 100 and the bond density p varies from 0.25 to 0.29. 
Curves are normalized by their values at p = 0.29. 



cally, due to the conservation of the energy in the system, 

(7) 



v 2 sr - 2 



where V and R are the external voltage and the bulk 
resistance of the system respectively, % is the value of 
the current on a bond, n, = 1 is the resistance of a bond, 
and the sum is taken over all bonds with nonzero current. 
For a unit external voltage source, Equation ([7]) can be 
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FIG. 6. (Color Online) Limiting current distributions 
f°°(i) = L _1 /_L(L _1 j) at different system sizes and bond den- 
sities. Colors represent the ratio L/(p—p c )~ o g as an indicator 
of the system size relative to the correlation length. 



rewritten as C = R^ 1 = Yl/b^bi wnere C is the bulk 
conductance of the network. Dividing by the number of 
current-carrying bonds N, we recover the conductance 
from the continuous current distribution: 



C 
N 



N 



i 2 fUi) di. 



(8) 



Equation (JSJ connects the scaling behavior of the bulk 
conductance to that of the current distribution. Since the 
second moment of the current distribution is dom- 

inated by the exponential tail of which is shown 

to be independent of the bond density p, the number 
of bonds carrying large currents (whose magnitudes are 
assumed to be larger than 10 -3 ) N{p) and the bulk con- 
ductance C(p) have the same scaling form with respect 
to the bond density p. The scaling of N(j>) with respect 
to p is shown in Figure El and the corresponding scaling 
form is: 



N(p) = (p - p c ) 



1.917±0.013 



(9) 



The balance of Equation ([5J conditioned on large cur- 
rents implies that there is an intrinsic consistency be- 
tween the power-law scaling of the conductance and the 
exponential tail of the current distribution. 

In conclusion, we have identified a robust and universal 
exponential distribution which persists above but close to 
threshold, describing the large current tail of the current 
distribution for the response of a cubic bond percolation 
network to an external field. In the supercritical regime 
above the percolation threshold, it is precisely this range 
of currents that is most relevant for describing and diag- 
nosing the macroscopic electrical response for materials 
applications. 
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